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Abstract. The finite- frequency dynamics of the 4 He superfluid phase transition can be 
formulated in terms of the response of thermally excited vortex loops to an oscillating flow 
field. The key parameter is the Hausdorff fractal dimension du of the loops, which affects the 
dynamics because the frictional force on a loop is proportional to the total perimeter P of the 
loop, which varies as P ~ a dfI where a is the loop diameter. Solving the 3D Fokker-Planck 
equation for the loop response at frequency u) yields a superfluid density which varies at T\ 
as kjVC'te- 1 ) _ This power-law variation with lu agrees with the scaling form found by Fisher, 
Fisher, and Huse, since the dynamic exponent z is identified as z — dn — 1. Flory scaling for the 
self-avoiding loops gives a fractal dimension in terms of the space dimension d as du = (d + 2)/2, 
yielding z = d/2 — 3/2 for d — 3, in complete agreement with dynamic scaling. 



1. Introduction 

The vortex-loop theory of the superfluid phase transition provides a very simple and physical 
picture of the transitionjTJElE]- An applied superflow acts to orient the thermally excited vortex 
loops, and their backflow currents screen the net superflow, reducing the superfluid density. As 
the temperature is raised towards the lambda point, larger and larger loops can be excited due 
to screening effects of the smaller loops, and for an infinite system at T\ a loop of infinite size 
drives the superfluid density to zero. This is not only the simplest renormalization-group theory 
of the transition, it also now appears to be the most accurate: two recent simulations [1] are 
in agreement with the loop prediction for the superfluid exponent v = 0.671688..., with the 
simulations finding v = 0.6717(1). 

This simplicity also carries over to the dynamics of the transition, where the characteristic 
phenomenon of "critical slowing down" can also be easily characterized in terms of the vortex- 
loop dynamics [5j. At low temperatures only small loops are excited, and with little frictional 
damping on them they are able to respond to a time-dependent flow field. As the temperature 
gets close to the transition, however, the larger loops being excited are subject to much stronger 
frictional forces, and can no longer keep up with changes in the applied flow, leading to a 
rapidly diverging superfluid relaxation time. The key parameter in this slowing of the dynamic 
response is the loop Hausdorff fractal dimension du, since the friction is proportional to the 
total perimeter P of a loop, which varies as P ~ a dtl where a is the average loop diameter. 

For the case of an oscillating applied flow at angular frequency co, the frictional force gives rise 
to a frequency-dependent diffusion length, where loops of diameter considerably smaller than 
this length can screen the applied flow, while the larger loops cannot. At the critical temperature 



the largest loops no longer respond to the oscillating flow, giving rise to a non-zero superfluid 
density at and above T\. This is a essentially a finite-size broadening of the transition, with the 
diffusion length now playing the role of the effective box size. 

The finite-frequency dynamics at T\ have also been formulated from dynamic scaling theory 
by Fisher, Fisher, and Huse [6]. We show here that vortex-loop dynamics give precisely the 
same results, and a dynamic exponent in full agreement with dynamic scaling. 



2. Vortex loop dynamics 

In the loop theory the renormalized finite- frequency superfluid density is given in terms 

of the zero-frequency density p s (0) by linear-response theory as [5] 
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where ag is the core diameter, the smallest loop size, and the response function g(ui, a) is found 
by solving the 3D Fokker-Planck equation [7] for vortex rings of diameter a, 
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where terms with derivatives of g have been neglected. Here D' is to within a constant the 
vortex-ring diffusion coefficient given by Donnelly [8] for a ring of diameter an. E is given by 

E = ir 2 K r — (ln(K' e ) + l) (3) 

where K r = h 2 p s ao/m 2 kBT is the dimensionless renormalized superfluid density, K = K r a/ao, 
and Flory scaling [3] gives!] 9 = d/((d + 2)(d - 2)) = 0.6 in d = 3 dimensions. At the critical 
point T = T\ an asymptotically exact solution of the loop recursion relations [H EJ [3] gives 
K r = K*ao/a where K* = 0.3875 is the fixed-point value of K. Since numerical evaluations 
show that this 1/a dependence is already valid at scale lengths of a few ao, the derivative in 
Eq. 1 is essentially a constant except for the smallest loops, and it is easy to see that E will 
then also be nearly constant except for the variation of K. However, at T\ K increases rapidly 
from its initial value Kqc = 0.309 at ao to its fixed point value K* = 0.3875 where it remains 
constant to the largest loop scales. Hence it is a very good approximation to take E « 5 to be 
a constant. In that case the integral of Eq. 1 is trivial, yielding for the frequency dependence of 
the superfluid density 

Rep s (w) -^/(^-i) (4) 

From further development of Eq. 1 it was shown in Ref . [5] that the superfluid relaxation time 
varies as 

/ a X^ -1 dH -i 

r ~UJ (5) 

where £ is the correlation length. From the definition of the dynamic exponent z, the loop 
prediction is then z = du — 1- Eq. 4 is then 

Rep s {u) ~ lu 1/z (6) 

1 The use in Ref. [3] of a 1/7? interaction between vortex segments means that their results are strictly only valid 
at d = 3. By using the more general interaction R~ < - d ~ 2 ' 1 the Flory minimization gives the vlaue of 6 used here. 
The fractal dimension however remains unchanged at Dh = (d + 2)/2. 



which is the main result of this paper, since it is exactly the d = 3 scaling prediction of Fisher, 
Fisher, and Huse [6j. The imaginary part, the dissipation, also has a universal dependence, 
which we express in terms of the phase angle 

\Rep s (oj)J 2 

which is also the same as in Ref.[6]. The value of z can be found from the Flory-scaling result 
|3j for the fractal dimension, dn = (d + 2)/2, giving then z = du — 1 = d/2 = 3/2, showing that 
the loop dynamics is in complete agreement with dynamic scaling [9j . Numerical simulations of 
the XY model have now verified this value of z to high accuracy [10J, finding z = 1.500(9). 

To check the approximations made in evaluating Eqs. 1 and 2 above, we have numerically 
integrated these equations, shown in Figs. 1 and 2. Only at very low frequencies is there a very 
slight difference from Eqs. 6 and 7, since at low frequencies a wider range of loop sizes is involved 
before g becomes small. 
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Figure 1. Full solution of Re p s (uj) 
from Eq. 1 (solid line) compared 
with a; 2 / 3 , Eq. 6 (dashed line). 



Figure 2. Full solution for the 
phase angle (solid line), compared 
with (p = 7r/3 from Eq. 7 



3. Discussion 

The agreement with dynamic scaling provides evidence that the Flory scaling result for the 
loop fractal dimension djj = (d + 2)/2 = 2.5 is quite accurate. Previous estimates of the 
fractal dimension from XY model simulations have found considerably lower values, however, of 
2.287(4) and 2.1693(4). These authors do not discuss the discrepancy between the two values, 
which is much greater than the claimed errors. The fractal dimension is extracted from the 
computed size distribution of the loops, a procedure which has been noted to be problematic by 
Olsson [12] . The problem is that on a lattice two loops that approach each other closer than a 
lattice spacing (which occurs with some frequency since the density is high near the transition 
temperature) cannot be resolved, so that there could either be two small loops or one large loop. 
A 50-50 guess is made in such a case, with the hope that the effect on the size distribution will 
average to zero. The simulations by Olsson gave evidence that this may not be the case, that 
the guess does lead to a distortion of the distribution. We note that a current simulation in 
progress not subject to this problem (the "worm" algorithm) is apparently now finding a value 



for the fractal dimension that is within a few percent of the Flory value of 2.5 (N. Prokof'ev, 
private communication) . 

Acknowledgments 

This work is supported by the US National Science Foundation, Grant DMR 05-48521. I thank 
S. Shenoy, N. Prokof'ev, A. Sudbo, and P. Goldbart for useful discussions. 

References 

[1] Williams G A 1987 Phys. Rev. Lett. 59 1926 

Williams G A 1999 Phys. Rev. Lett. 82 1201 

Williams G A 2004 Phys. Rev. Lett. 92 197003 
[2] Shenoy, S R 1989 Phys. Rev. B 40 5056 

[3] Chattopadhyay B, Mahato M, Shenoy S R 1993 Phys. Rev. B 47 1515 

[4] Campostrini M Hasenbusch M Pelissetto A and Vicari E 2006 Phys. Rev. B 74 144506 

Burovski E Machta J Prokof'ev N and Svistunov B 2006 Phys. Rev. B 74 132502 
[5] Williams G A 1993 Phys. Rev. Lett. 71 392 

Williams G A 1993 J. Low Temp. Phys. 93 1079 
[6] Fisher D S Fisher M P A and Huse D 1991 Phys. Rev. B 43 130 
[7] Donnelly R J and Roberts P H 1971 Phil. Trans. Roy. Soc. 271 41 

[8] Donnelly R J 1991 Quantized Vortices in Hell (Cambridge University Press: Cambridge) 

[9] Hohenberg P C and Halperin B I 1977 Rev. Mod. Phys. 49 435 
[10] Oh S K, Yoon C N, Chung J S and Kang H J 2004 J. Korean Phys. Soc. 44 491. 
[11] Hove J, Mo S and Sudbo A 2000 Phys. Rev. Lett.85 2368 

Camarda M, Siringo F, Pucci R, Sudbo A, Hove J 2006 Phys. Rev. B 74 104507 
[12] Olsson, P 2002 Europhys. Lett. 58 705 (Preprint |c"ond-mat/0103345[ ) 



